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Abstract 

In this paper, we investigate the performance of selection cooperation in the presence of imperfect 
channel estimation. In particular, we consider a cooperative scenario with multiple relays and amplify- 
and-forward protocol over frequency flat fading channels. In the selection scheme, only the "best" relay 
which maximizes the effective signal-to-noise ratio (SNR) at the receiver end is selected. We present 
lower and upper bounds on the effective SNR and derive closed-form expressions for the average symbol 
error rate (ASER), outage probability and average capacity per bandwidth of the received signal in the 
presence of channel estimation errors. A simulation study is presented to corroborate the analytical 
results and to demonstrate the performance of relay selection with imperfect channel estimation. 
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I. Introduction 

The promise of high spectral efficiency and capability of providing great capacity improve- 
ments in a wireless fading environment, as reported by [DL1 and 0, has led to widespread 
interest in multi-input multi-output (MIMO) communications. However, due to size, cost, and/or 
hardware limitations, a wireless device may not be able to support multiple transmit antennas. 
Cooperative diversity was proposed as an alternative to MIMO systems. It has been demonstrated 
that cooperative diversity provides an effective way of improving spectral and power efficiency 
of the wireless networks O, H. The main idea behind cooperative diversity is that in a wireless 
environment, the signal transmitted by the source (S) is overheard by other nodes, which can 
be defined as "relays" [5]. The source and its partners can then jointly process and transmit 
their information, thereby creating a "virtual antenna array", although each of them is equipped 
with only one antenna. It is shown in ll6l-[fT4l| that cooperative diversity networks can achieve 
a diversity order equal to the number of paths between the source and the destination, however, 
the need of transmitting the symbols in a time division multiplexing (TDMA) fashion limits the 
improvement in capacity. Additionally due to the power allocation constraints, multiple relay 
deployment is not economic. Relay selection then was proposed to alleviate the loss in spectral 
efficiency caused by multiple relay schemes and also to moderate the power allocation constraints. 

Recently, there have been considerable research efforts on the performance analysis of coop- 
erative diversity including the derivation of closed-form formulas for the average symbol error 
rate (ASER) 0, [H, [H, El-flE!, 031, the outage probability 1171 191 ifTOl, ifBl, |fl6l, and 
the average capacity (ergodic capacity) flTJ, IfTOl . |fTfl. ifTTl . Most of the current works, however, 
assume the availability of perfect channel state information (CSI) at the destination (D) terminal. 
Channel estimation for relay communication has also been studied in literature as in lfl8l . [fl9l . 

Related work and contributions: The impact of channel estimation error on the ASER perfor- 
mance of distributed space time block coded (DSTBC) systems has been investigated in ll20ll . 
assuming the amplify-and-forward protocol. Building upon a similar set-up, Gedik and Uysal GTI 
have extended the work of [20] to a system with M relays, assuming the amplify-and-forward 
protocol. In 11221 . the symbol error analysis is investigated for the same scenario as in 11211 . In 
this paper, we investigate the effect of channel estimation error on the ASER, outage probability 
and capacity of a cooperative diversity system with relay selection. We first introduce our system 
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model in Sec. HO In Sec. UTTl we explain the Gaussian error model introduced by ||23l and use it 
in the selection scenario to obtain the maximum ratio combiner (MRC) output SNR. In Sec. |V] 
we study the ASER performance of our system in the case of imperfect channel estimation. We 
develop upper and lower bound analysis for the general relay case to formulate the relay selection 
scheme. We then derive a closed-form formula for the ASER performance with imperfect channel 
estimation. We also derive closed-form formula for the ASER performance of our system in the 
high SNR regime as introduced in 0, [|24]|. We also study the outage probability and average 
capacity per bandwidth in Sec. |VI]and Sec. IVIIl respectively. Diversity analysis is also provided 
in IVIIIl Simulation results are presented in Sec. [DD followed by conclusions in Sec. |X] 

Notation: E{-}, {■}* and x, denote the expected value, the complex conjugate and the estimate 
of the variable x, respectively. 7 t lb and j i ub stand for the i th lower bound and upper bound 7, 
respectively. 7i* indicates that the relay i with SNR 7 is selected. f x (x) and F x (x) denote the 
probability density function (PDF) and the cumulative density function (CDF), respectively. 



In this section we consider a system in which a source node S transmits information to a 
destination node D with the help of the best relay node which is selected under a defined 
criterion from M available relays. The transmissions are orthogonal, either through time or 
frequency division. 

In the first data sharing time slot, the source node communicates with the destination as well 
as the relay nodes. In this phase, the signals received by the destination and each relay are 



where x, y sd , and y si denote the transmitted signal with unit energy and the signals received 
at the destination and the i th relay node, respectively. h si and h sd are the channel coefficients 
of the source-relay and the source-destination channels, which include the effect of fading. £ sd 
and £ si represent the average signal energies received at destination and the relay i terminals, 
respectively, taking into account the path loss and shadowing effects. n sd and n si are additive 
white Gaussian noises (AWGN) in the corresponding channels with the same variance N , i.e., 



II. SYSTEM MODEL 




(1) 




(2) 



n 



'ad' 



n 



CAf(0 : N o ). 
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In the next time slot each relay node normalizes its received signal and retransmits it to 
the destination. For the i th relay, the normalization factor is \j £ sl \h ai | 2 + N and the signal 
transmitted from the this relay is Q 



\/S h x + n . 

V SI SI ' St /' m >\ 

Xi = ■ (3) 



'S.\h ] 2 + N 

SI \ SI \ u 

Since the relays practically do not have a perfect knowledge about the channel h si , i — 1, . . . , M, 
the power normalization procedure in each relay is performed based on the estimate of the 
corresponding channel. We assume in this paper that E{|x| 2 } = 1. Based on the signal 
received by the destination from the i th relay node is Q 

= \ZS~ d h. d Xi + n. d , 



JsjT d 

KiK x + n ui ( 4 ) 



= aih h,x + n. ,. (5) 

where h. d is the channel gain from this node to the destination. £. d is the average signal energy 
received at the destination via the i th relay in its corresponding time slot which accounts for the 
shadowing and path loss effects and oa = v ld . n. d ~ CAf(0, N ) denotes the AWGN 

V £ »\ h si\ 2 + N 

of the relay destination channel. n id is the equivalent noise term in y. d . It can easily be shown 
that, conditioned on the channel realizations n. d ~ CJ\f(0,ufN o ), where 

2 = 1 + *d^d\ _ (6) 

£\h\ 2 + N 

S% \ SI I 1 u 

Supposing only the i th relay participates in cooperation we look for the condition under which 
we can select the best relay, Ri*, by searching RiS for i = 1, . . . , M. 

Since each relay only amplifies the signal from the source, the destination is the only place 
where an estimate of the information symbol x is computed. In order to achieve maximum 
likelihood performance, the signals from both diversity branches (direct S — >■ D branch and 
the branch via the i th relay) are combined using a maximum ratio combiner (MRC). Decoding 
of x is delayed until the relayed signal containing the information symbol x is received at the 
destination. Since the noise power is not the same on the two sub channels, both diversity 
branches must be weighted by their respective complex fading gain over total noise power on 
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that particular branch before the combiner. Thus we obtain the estimated information symbol as 

h*,/£T a*h*.h* 



where 



SZ ! 


h . 

SI 


2 


1 %d 1 


h. d \ 

id 1 


2 + N 




h .1 


2 


+ N 



/\ si \ si \ id \ id \ u ZO\ 

"" = V £.|&.|' + «. ' W 

V si I sx \ u 

since the destination only uses the estimate of the relay-destination channel. 

III. Modeling of channel estimation error 

Let the true channel gain and its estimate be h ~ CAf(Q, crfy and h ~ CjV(0, ct|), respectively. 
It is shown in ll23l . Il25ll that it is possible to consider the channel estimate as 

h = h + e, (9) 

where e ~ CjV(0, erf ) is a zero mean complex Gaussian noise. As h and h, are jointly Gaussian, 
the conditional probability density function (PDF), f(h\h), would also be Gaussian with the 
mean B26l 

2 

E{h\h} = \h, (10) 

and variance of 

2 

var{%} = (1 - (11) 

By (flOl) and (TTTb and assuming a least squares (LS) estimator, we can consider the channel as 

h = ph + d (12) 

where p = (a%/cr?) = P h h{°hl 'h) Phh * s me correlation coefficient between h and h and 
d ~ CjV(0, erf,). In this paper, erf, d , crfj si , crf, id , p sdl p si and p id are the channel estimation error 
variances and their corresponding p factors of the S — D, S — Ri, and R4 — D links, respectively. 
According to (fTT|) 

<?d = 0--p)°h 

= (P~P 2 K (13) 
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Using (PT21) . we can rewrite the first term of © as 



D ° d = Sd N ^ [\^(Ps d h s d + d J x + n s 



x + ±_g_«(, £d + I (14) 



10 

Here, we have inserted (fl"2l) in (OQ) then the resultant expression is substituted into ©. Following 



[|22|. the received direct path signal after the MRC can be decomposed into the message part 
given as Psd£ ^ sd and the noise part which is given as y/ ^^ sd ( \f£~ d d sd + n sd ). Finally, due 
to the independence of the h sd , d sd and n sd , the effective SNR due to the direct path is given by 

E x {|E{D d }| 2 } 



7 



,s- d 



™{D sd } 

tP 2 Jk d \ 2 pI% 



1 + £sd a 2 1 + e , ' 

1 ' N ° D,sd ad 

where 7 sd is the estimated received SNR due to the direct path and e ad = jy^crf, ^. 
Following the same approach, we can write the second term of © as 



(15) 



ctih* h* 

si id 



^i{p sl h sl + dj{p. d h. d + djx + n. o 



(16) 



ufN 

We can extract the message part, error part due to channel estimation error, and the noise part 
of the received signal a|j] 

a 2 

M = -^r—p P-Jh .\ 2 \h.,\ 2 x, (17) 

,",2 AT ' S1 ' zd \ st \ \ zd \ ' V * 

LO i iV 
2 

.\*h*d.,+p. „\h. fh*d . 

A 2 AT si I si I zd 1 <id I id I S J si 

iV 

+h* h* d d jx, (18) 

^=-hV 4 n 4(J , (19) 

respectively. To obtain the effective SNR expression, we need to find the ratio of the signal power 
to the overall noise power. Alternatively, the received SNR can be calculated as E "^jg^j |2} 11271. 

'Conditioned on h gi and h id the distribution of d 1 8 b is approximated as Gaussian. The validity of this approximation has 
been confirmed by simulation. It should be further noted that the same conclusion has been reported by 1221 and the references 
therein. 
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Since h ., h ,, d and d., are zero mean independent processes, we can write the effective SNR 
of the selected path as 

£ o 2 \h \ 2 S o 2 \h I 2 

_,eff at' si' si\ idr id \ id\ (20) 

i /V" 2 C£2i|/j |2f1 4- n 2 -^n 2 ) 4- -«M/j l 2 (1 4- n 2 ^rr 2 ) 4 ^-rr 2 ^-rr 2 4-11" 



With a simple manipulation, we may write (1201) as follows 

„,eff ^ai^ id^silid 1 \ 

where A . = 1 + p 2 e.,, A.. = 1 + p 2 e ., 7 . = £s '^ s ' ;l , 7 , = ^|M-, e . = £st ° f f ,ai , and 

r si id> id r 1(i si' Isi No 1-1 ' ld iVo 'si No 

e id = id N ° M ■ Following J22l . 1 + e si e id can be ignored □; hence, (|22l) can be approximated as 



U, 122 



eff _ P s iP id'Ysiiid 



^ett ' si 1 id 1 al '^d (22) 

i 7 .A . + 7. .A.j 

I sz si 1 /id %d 

IV. selection Strategy 

In the selection based amplify-and-forward scheme studied in JU-JU, J28l . [|29ll the best relay 
selected by the destination terminal is the one that leads to the maximum total received SNR, 
7 r , which is given by 

7r = € + 7i?j (23) 

where 

7 f = arg max{ 7 f }, (24) 

for i = 1 . . . M. 7 r , which is the total received SNR at the destination, is the addition of the 
SNR received from the direct path, 7^, and the SNR received via the selected relay path 7??. 
The selection steps are mainly as follows. The source terminal sends a ready to send (RTS) 
packet to all the relays and the destination. As the relays receive this packet they estimate their 



corresponding source-to-relay channel h s i. They divide their signal by the factor y 8 S i\h si \ 2 + N 
and wait for the clear to send (CTS) packet from the receiver side. As each relay receives the 

2 For a fixed SNR, i.e., ^ = 5dB, we have noticed that the mean square error MSE between ji S in fT\\ and f22\ is 10 10 

off 2 

percent of y\ for y si — y id — 30dB. 
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CTS overhead, it estimates its corresponding relay-to-destination h id channel and will start a 
timer which is given by 

1 

where 7 2 ? ff is given by (|2D) . Each timer T{ counts down to zero right after receiving the CTS 
packet. The first timer that reaches zero sends a flag to all the relays and the destination, 
introducing itself as the best relay. Then all the relays would remain silent, and the selected 
relay forwards its signal to the destination, with an overhead package containing information 
about its source-to-relay link h si and the relay-to-destination channel h^. It is also possible for 
the relays to use a separate log 2 (|/isi| 2 ) + log 2 (|/ii<i| 2 ) bit feedback link to send the CSI to the 
destination. Further details on the selection procedures can be found in lfl4l . lfT6l . 

V. ASER ANALYSIS 

In this section, we derive an ASER expression for the relay selection scheme. We first introduce 
the following theorems 11301 

• Theorem T. If X and Y are two random variables with Y = ml, then 

M7) = Afx(-). 

\m\ m 

• Theorem 2: For two independent random variables X and Y if Z = min(X, Y), then 

fz{l) = fx{l) + fv{l) - fx{l)F Y {l) - fy{l)F x {l). 

• Theorem 3: Let Y = maXj(Xj), i = 1, • • • M, where XjS are independent random variables, 
then 

M 

Fy( 1 ) = l[F x Xl). 
i=l 

It is shown in 0, [H, d, [H, d-El ESI that the symbol error rate (SER) in a cooperative 
scenario is given by Q(^/k%), where 



1 r°° -ft 
Q(x) = —= I exp(-— ) dt, 



/27T J x 2 

and k depends on the modulation schemed 

3 For instance for BPSK modulation k = 2. The analysis can be generalized to QPSK or QPAM modulation as well [5j. 
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To find the average symbol error rate, we need to integrate Q(-^k%) over the probability 
density function (PDF) of %, i.e., 

POO 

P e = / Q(y/H) d%. (25) 



Jo 

Finding the exact expression for / 7 can be quite cumbersome, therefore, in the following and 
to simplify analysis, we develop lower and upper bounds on %. It is straightforward to show 
that for two arbitrary independent random variables, the following inequality holds [|6l 

1 XY 

- min (X,Y)< -j^-y < min (X, Y). (26) 
Applying d26j to <yf in d22j) yields 

7 M6 < if < 7^, (27) 

where 

P 2 P 2 

li,u b = j^f- min (%Av 7.AJ, (28) 

si id 

7 M6 = ^^min(7 si A si)Ti AJ> (29) 

si id 

Finally, j r , can be bounded as 

l lb <lr<l^ (30) 

where 7 ;b and j ub are the lower bound and the upper bound SNR values defined as 

p 2 ff Bd 

7 U (, ~~ i T I" 7, ; *„ b ; (31) 
7, = :f2T L + 7 4 *„, ( 32 ) 

respectively, where 

%*u b = max {7 l( ,J, (33) 
7i* B = max {7 M J- (34) 
Using (HI, CD) and ([32]), (f25]) can be bounded as 

P lb <Pe< P ub , (35) 

where P lb and P ub , are the lower and the upper bounds for ASER's. In this paper, we limit our 
discussion to the lower bound analysis. The upper bound analysis can be obtained similarly. 
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Since the ASER is a decreasing function of SNR, a lower bound for J °° Q{y/krf r ) fy r d% would 
simply be J ° c ' Q(y/kl7 b )f^ b dnf ub . 

Since our final goal is to obtain the PDF of 7 . we need to find the PDF of 7. ,, for i = 

O i ub • 1 ,ub 

1, . . . M, first. Note that the PDF of 7 sj can be expressed in terms of the average SNR j si = 
ai }j 3i as fxf (7) = ^-exp (^). Similar expressions can be found for the PDF of 7 sd and 

s * ^si ' 'si [""I 

7 id as well. Thus, using Theorem 1, the PDF of r ) iub can be written as □ 



A,. l W=^exp(-^), 



(36) 



where = ^ SJ = id ^ si ^ d • Using Theorem 2 and noting that 7.^ 



A .7 .+A. ,7. j ' 



maxj (7 lu J, we obtain 



c?7 

«n£i^..(7) 



c?7 



A/ 



* = 1 



The closed-form formulation for / 7 ^ (7) would then be (see Appendix I) 



(37) 
(38) 
(39) 



M \PJ 

A„ b (7) = ££(-i; 

p=l m=l 



exp(^)-exp(-7[^b] m ) 



where // = ^ and \I/ P is a binary permutation matrix with dimension of J x M which is 
defined in the Appendix I and 



(40) 



b^ii...^ 



L /?i 02 0. 



M 



(41) 



Having noted that 



P„ 



0(7) = i-iErf(^). 



5 " §M^)] 



(42) 



(43) 



4 Here we have to note that ■y i u6 s are not i.i.d random variables, since /3jS are not analogous for different relay paths. 
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by substituting (|40l ) into (|43l) and setting j ub = 7 2 , we can write P lb as 



fc/2 



[^Pb] m +fc/2 



(44) 



Here, we have used the fact OT 



/ 7 Erf(a 7 )e X p(-6 2 7 2 )d7 = — -== b^O. (45) 
Jo 2b 2 yb 1 + a 2 

Note that the upper bound analysis can be derived similarly. 

To gain further insight into the performance of the selection scheme under consideration, we 
resort to the high SNR analysis. Defining 7 ub = r y ub / , y ub , where j ub is the average SNR, and 
following Q, 11241 . we can approximate the asymptotic behavior of /* in the high SNR regime 
with its Mclaurin series expansion, i.e., 



f% b = a Z + o(\ b ), (46) 
where a = ~ d Z ub (0), if the derivatives of up to order n — 1 are null and 

ub 

t=n+l "ub 

It can be easily shown that by replacing (l46l) into (1431) . the asymptotic behavior of ASER 
(7^ -» 00) is given by ED 

p _ n" m t 1 i(2m-i ) 1 g ! ^ 

ifc 2(n + l)A;(«+ 1 ) n! <9 7 ™ 1 J ' V ; 

where we have used the fact that f~ (7 ,) = ^-/^ (^M. 

Similarly, the asymptotic upper bound ASER can be found as 

p , n^m-l) ! ^/, 

2(n + 1)^ 1 ) n! <9 7 ™ 1 J ' 1 ; 

Theorem 4: For all n < M, the values of 9 ^ b (0) and g ^ (0) are zero and for n = M 

ub lb 



P, -> 



nS(2m-l) 1 



2(M+l)M M+1 )/i^/3.. 



' (M + l)£;(^+i) i/ll A 
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d 'ad 



where ^ = j 

Proof of Theorem 4: (see Appendix II) 

VI. Outage probability Analysis 

The outage probability P out , which is a valid measure of performance in slowly fading 
channels, is defined as the probability that the source-destination mutual information X sd falls 
below the transmission rate R. Defining 

Z d = \\og 2 {l + lf d +lf)<R, (5D 
the outage probability is given by 

Pout = P(|l0g 2 (l+ 7r ) <R) 

= P(%<2 2R -l). (52) 

Substituting the introduced upper bound for 7 r leads to adopt an upper bound on the outage 
probability, i.e., 

Pout<P^ t = P(l ub <2™-l) 

= ^(2^-1). (53) 



We can simply find F 7 by integrating (|40l) with respect to 7 the result of which is given by 



7 U 

( M 

M \P 

P-rJ-T) = E E(" 1 ) (P+1) [l + exp(- 7 [^b] m ) - JSfc T exp(^ ) 

p=l m=l 

(54) 

We can apply the same procedure for the lower bound analysis. 

VII. Average Capacity 

The maximum achievable transmission rate under which the occurred errors can be recovered 
is called the system ergodic capacity. The average capacity per time slot for S-AF scheme for 
the perfect CSI case is defined as ll32l 

C = T So lo&(1 + %)f% {lt)dlt (62) 
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in bits per time slot, where W is the transmitted signal bandwidth, and j t is the total received 
SNR assuming perfect CSI. To obtain a formula for the average capacity with imperfect CSI, 
we have to consider an information theoretic approach. Regarding (fl4l) and (fT6b . the total signal 
in the MRC output at the destination node after two time slots is 



y o = h Q x + n , 



where y o = 
scalar, and n r 



sd i* d\ 



rsd sd I sd I 



sd s d 
Nn 



N 

EA^x + n. 



M 



(63) 



x is the transmitted signal and is a 
, where V and M are given in (TT8T) 



and (fT9l) , respectively. Normalizing the received vector y G by the noise variance we arrive at 

y = ha; + n, (64) 



where 



where 



sd 


D i*d 






P sd ^ sd ' l sd 






M 







l T 



(65) 



(66) 



and n ~ CJ\f(0, 1) . If the noise vector n is independent of the input signal x, then the capacity 
of the system in (|64l ) is given by 



C = arg max X(y, x) 

= arg max H(y) — 1-L{y\x) 
= arg max H(y) — T-L(n\x) 

= arg max H(y) — K(n) 

/(*) 

W 



log 2 (1 + h^h) . 



(67) 
(68) 

(69) 



where %(•) denotes the differential entropy and f(x) is the PDF of the input signal x. In our 
system the noise vector n depends on the input signal x, and on the other hand, the noise vector 
n is not Gaussian. Therefore, (|67|) will not lead to (1681) and hence we are not allowed to use the 
capacity formula in (l69l . Instead, using the following inequality 11331 



H{n\x) < H(n) 
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we can bound the capacity with the worst effect of the additive noise on the capacity of the 
system [|34l . which is when the noise vector n is Gaussian. It has been shown in [|35l . Theorem 
1, that if the transmitted signal and additive noise are uncorrelated, the worst case noise has 
zero-mean Gaussian distribution with the same power of the additive noise. Furthermore, it can 
be easily seen from (|64|) that E{nx\h sd , h si , h id } = 0. Therefore, the additive noise n can be 
replaced by a Gaussian noise with the power 1, leading to the worst case capacity given by (|69l) . 
By substituting D sd and D. d from ([14)) and (fT6l) in (l66l) and constituting (l69l) . we have 

h^h = 7r . 

Therefore, the worst case average capacity is given by 

W f°° 

C WO rst = — / log 2 (l +7r)/ 7 r(7r)^7r- 

n 

To obtain the lower bound average capacity^ C lb < C W0TSt < C, we substitute 7 r with j ub and 
f 7r with (|40l) . Therefore, we get 

M 



W 

p=l m=l 



(70) 



Here, we have used the fact that [[3111 

f°° 1.44 

/ exp(— cry) log 2 (l + 7)g?7 = exp(a;)-E , z(— a) 

Jo " a ' 



1.44 

= exp(a)El(a), 

a 

where £'1(7) = —Ei(-'y), and 

EiOy) = / — exp(— a;)<ir 
J 7 a; 

is the exponential integral function. 

5 log 2 (l + 7)/7(7) is an increasing function of 7; hence, for j ub > y r , / °° log 2 (l + J ui )f~f ul (7 ui )dj ul < J °° log 2 (l + 
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VIII. Asymptotic Order of Diversity 

We shall next analyze the underlying selection cooperation system from the diversity point of 
view. The diversity order is given by the magnitude of the slope of the outage probability as a 
function of SNR (on a log-log scale) [4j. Equivalently a system achieving a diversity order d(R), 
at transmission rate R, has an error probability that behaves as P ub (SNR) oc SNR~ d ( R ^ at high 
SNR [4J. Utilizing the high SNR analysis which is indicated in ([50]), the ASER is implicitly 
proportional to 7 sd ,7 s4 , and 7 id . Assuming 7 sd = «i7 si = a 2 7 ld = P , then we have 



p 2 p 11 p2 p 2 p 2 



p 2 P Q J- J- P p 2 p 2 

"sd i— i or al r id 

t!1 — ^ — • ( 64 ) 



oc 



oc 



p 2 pM+1 11 p 2 p 2 

This implies that the diversity order in a selection cooperation scenario is highly affected by the 
channel estimation error. In the case that the channel estimation error is zero or it reduces by 
increasing the received SNR, e sd , A 8i , A id , p sd , p si , and p. d , for % — 1, . . . , M, would all be equal 
to unity at high SNR regimes; hence, 

P«{Po)«(h U+l - (65) 

Thus, the maximum diversity order is achieved in selection cooperation schemes with no channel 
estimation error or in the cases that the channel estimation error reduces with increasing SNR. 
On the other hand, if the channel estimation error is independent of the received SNR, noting 

that e sd = ^°D,«l> 6 * = No a D,si> ^ = l^ a2 D,id> X s, = 1 + Pi ^ and \ d = 1 + Pl e si> We Can 

replace e sd and similarly A si , and A id by P , at high SNR; hence, from (l64l) the diversity order 
of the system with selection cooperation tends to zero as P , approaches infinity. This means 
that in such a system we expect an error floor in the ASER or outage probability curves at high 
SNR, in which the performance of the system would not improve by increasing the SNR. 

IX. Simulation Results 

In this section, we consider two relay selection scenarios with M = 2, and M = 3, respectively. 
The transmitted symbols are drawn from the PSK constellation where we assume k = 2 and 
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the received energies are assumed to be (£ sd ,£ si ,£ id ) = (1,0.5,0.5) x P for all relays. It is 
assumed that the power budget of each relay can not exceed Pq/2 and once a relay is selected, it 
operates with full power and the other relays remain silent. All the channels are i.i.d zero mean 
and unit variance complex Gaussian random variables. For convenience, we assume the special 
case, where the channels estimation errors of all links have the same variance, i.e., a 2 D sd = 
a D,si = a Did = a D an d all the correlation coefficients are equal, meaning p sd = p s . = p. d = p. 
The variance of noise components is set to N = 1 and R = 1 bps/Hz. Here, unless stated, we 
assume that a 2 D is independent of the received SNR. 

ASER analysis: In Figs. [2] and [3l we confirm the validity of the derived expressions for the 
cases M — 2, and M — 3, respectively. As the figures show, the derived lower bound for the 
ASER is fairly tight. The high SNR analysis also converges to the analytically derived lower 
bound ASER curves at high SNR. As it is obvious in these figures, if the channel is estimated 
with even a very small error, an error floor is observed at high SNR. Increasing p yields a better 
performance in the ASER of the system. Also we can deduce from these figures that selection 
among a larger number of relays improves the ASER performance of the system. 

Outage probability analysis: Figs. |4] and [5] considers the credibility of our analytic results for 
M = 2, and M = 3. In these figures the derived formula for the outage probability is compared 
with the simulation results. A fair match is observed and the same error floor is seen for the 
case that the system faces a non-decreasing channel estimation error with the received SNR. It is 
observed in Figs. |4] and [5] that the selection scheme with larger number of relays outperforms the 
scenario with less number of relays in terms of outage probability. Also the effect of increasing 
of p and thus, decreasing of a 2 D in the improvement of the system outage probability is shown 
in these figures. We further study the effect of relay location in selection cooperation. In this 
scenario the distance between S and D is normalized to one. The received energies are modeled 
as £ . — and £ ., = ° ,5J ?° where d is the normalized distance between S and FL. 

Fig. [6] shows the outage probability of a relay selection scenario with M = 2, versus d. The 
two relays are located between the source and destination nodes with the same distance of d 
from the source node. As it is noticed the lower bound analysis curves fairly match the curves 
obtained from simulation. The lowest outage probability is achieved when the selected relay is 
in the middle of the line of sight of source and destination. 

Asymptotic Diversity Order: Figs. [2H2] also considers the credibility of our diversity analysis 
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for M = 2 and M = 3. The slope of the curves associated with p = 1 in these figures imply the 
diversity orders achieved by the selection cooperation scenario. Diversity order M + 1 is seen 
for the perfect channel estimation scenario. Occurrence of error floor in the curves associated 
with the imperfect channel estimation accounts for the diversity order of zero at high SNR's 
which is analytically obtained. 

Fig. |7] shows the ASER performance curves for S-AF and all participate amplify-and-forward 
cooperation (AP-AF), where all the relays participate in cooperation. It is noticed that both 
scenarios exhibit the same slope, i.e., same diversity order. This means that S-AF as well as 
AP-AF in the cooperative system achieve full diversity order, i.e, M + l, where M is the number 
of relays. Although selection degrades the ASER performance of the system, it improves the 
throughput of the system significantly. Also S-AF cooperation requires less power consumption 
than the AP-AF cooperation scheme. 

The same results are achievable for the case that the estimation error variance is decreased by 
increasing the SNR. Fig. [8] shows the simulation results for M = 2 and M = 3 simultaneously. 
In this figure the received energies are assumed to be (£ sd , £ si , £. d ) = (0.5, 0.5, 0.5) x P , for all 
relays. Diversity orders of 3 and 4 are achieved for M = 2 and M = 3, respectively. 

Average Capacity: Fig. [9] shows the capacity of the S-AF scenario when we select the best relay 
out of M = 2, relays. A capacity ceiling is observed when the channel is estimated with error. 
After the effect of the receiver noise is completely diminished, which occurs at a certain SNR, 
the capacity improvement is insignificant from that point on. This is in fact because increasing 
the SNR does not reduce the channel estimation error in our scenario. As it can be seen, the 
capacity of the system is very sensitive to the channel estimation performance. A fair upper 
bound is observed in our simulations. Furthermore, Fig. [10] shows a comparison between S-AF 
and AP-AF cooperation for M = 2,3, 6, and 9. Obviously, S-AF cooperation exhibits a huge 
average capacity improvement specially for m > 2. Selection out of larger number of relays 
would slightly improve the throughput of the system. 

X. Conclusion 

In this paper we study the effect of channel estimation error in cooperative systems with relay 
selection. We analyze the ASER and outage probability performance the system in the presence 
of estimation error. We further derive closed form formulations for the average capacity of the 
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selection cooperation with imperfect channel estimates. Our analysis shows that the selection 
cooperation scheme still preserves the full diversity order. A simulation study is presented to 
corroborate the analytical results and to demonstrate the performance of relay selection with 
imperfect channel estimation. 

Appendix. I 

Deriving the Probability density function of 7.* . 
Consider a case with three relays, i.e., M — 3. From Theorem 3 we have 



^W = II( 1 -«P(-^))- 



Using d3U) to d3~9i we have 



W7) - E^ ex P(-J) fl (l"«P(-i)) 



1=1 



i = 1 



exp(-£) - exp(- 7 (£ + £)) - exp(- 7 (i + £)) - exp(- 7 (£ + ^ + ^)) 
exp(-£) - exp(- 7 (4- + 4-)) - exp(- 7 (4- + 4-)) - exp(- 7 (4- + 4- + 4-)) 



ft 
1 

+ s 

1 

+ A 

= ^ exp( -^ ) + ^ exp( -i ) + ^ exp( -S ) 



exp(-£) - exp(- 7 (i + i)) - exp(- 7 (i + £)) - exp(- 7 (i + i + £)) 



+ ( ^ + ^ + ^ )eXP( - 7( ^ + ^ + ^ )) - 



(66) 



Let \l/ p be a binary permutation matrix with dimension of 
possible p combination of M binary bits. In this example 





1 










1 


1 










1 





,^ 2 = 





1 


1 










1 




1 





1 



We also define b 



1 T 



x M. Each row shows one 



1 1 1 



and let [\I/ p b] TO be the m row of the vector \l/ p b. Using 
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the definitions we just pointed out we can write (1661 ) as 



A^(7) = E E(- 1 ) (P+1) t Vl/Pb W^P (-7[^ p b] m ) . (67) 

p=l m=l 

With the same procedure for M relays, using (l36l) and substituting F 7 fc = [1 — exp(— into 
(|39| ), we obtain 



= E E^ 1 )^ 1 ^]™ 6 ^ ("7[^b] m ) , (68) 

p=l m=l 

where b is defined in (I4TT) . 

Furthermore, it can be easily shown that the PDF of ji = ad is given by ^(7) = -4 exp(— 
where // = /f , and the PDF of 7 u6 , can be found as follows 



^(7) = ^(^<7) 

= + 7^ < 7) 



5 yy 

= fry Jo P ^- 1 ' 1 ^ )f ^ dl ^ 

= dllo ^ (7_7 ^ )/ ^>^ 7 ^- (69) 
Using ^(7) = 1 — exp(— V) and inserting ((681) in (l69lh we obtain 

A* (7) = E E(" 1 ) (P+ V / [1 " ex P (-^^)][^b] m ex P (- 7( ,J**b] m ) ^ 

p=l m=l " 



/Z[*Pb] 

p=l m=l 1 J 



(70) 



Appendix. II 
PROOF of Theorem 4 

From the basic principles of the moment generation function (MGF) of 1 ub = ^ + %* ub , we 
have 

M 7 (a) = M lt (8)M 7 (a), (71) 
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thus using the initial value theorem, the value of Q r y b (0) is given by 

ub 



Qn f 

— -^(0) = lim s n+1 M~ (s) 

lim sM (1 (s)s n M Ti (s) 



Qn-1 j 

^^(0). (72) 



where 

9 n "V 7 «9 n F 7 

(0) = (0) . (73) 

i*ub i,ub 

In light of Theorem 3, we can simply deduce that for n < M, 8 n * (0) is the summation 

i*ub 

of some terms in the product form, where each product form contains at least one cumulative 



density function (CDF) of ^f, at zero. Since 1^(0) = 0, we have }* ub (0) = 0; hence, 

' i*ub 

from (E3), -§^■(0) = 0. For n = M, we have 

ub 



d M ~ x h d M F^ 

i*ub i*ub 

M 

= MlUf^ + e^J, (74) 
i=i 

where in this equation e(7.* u J denotes all the other terms. Note that each term in e(7 jiu J has 
a product form and contains at least one CDF of 7^ 6 at zero and e(7.* J = 0. Substituting (1741) 
into (22) yields 

-q^(0) = MlU0)l[f^0) 

tub %=\ 

M 



(75) 



By (l28l) . (|29l) and Theorem 1 and following similar steps to above, we obtain 

QMf 1 M 

^#(0)^M!-n_. (76) 



Finally by substituting d75j) and <|76j) into (|47j) and (gSj), we obtain the result in <|49j), which 
completes the proof. 
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Fig. 2. ASER vs. SNR. Relay selection scenario, M=2. 
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Fig. 3. ASER vs. SNR. Relay selection scenario, M=3. 
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Fig. 4. Pout vs. SNR. Relay selection scenario, M=2. 
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Fig. 5. Pout vs. SNR. Relay selection scenario, M=3. 
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Fig. 6. Outage probability for M = 2 relays vs. distance for different SNR values. 



January 4, 2011 



IEEE TRANSACTIONS ON WIRELESS COMMUNICATION, SUBMITTED DEC. 2010. 



29 



10 




5 10 15 20 25 30 35 



Po(dB) 

Fig. 7. ASER vs SNR. Performance of S-AF and AP-AF Cooperation for M=2. 
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Fig. 8. ASER Versus SNR. S-AF for M=2 and 3, when the estimation error variance decreases with increasing the SNR. 
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Fig. 9. Average capacity per bandwidth vs. SNR. Relay selection Scenario, M=2. 
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Fig. 10. Average capacity per bandwidth vs. SNR. Comparison between S-AF and AP-AF cooperation schemes. 
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